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ABSTRACT 

We have identified five new R Coronae Borealis (RGB) stars in the Galactic bulge 
using the MAGHO Project photometry database, raising the total number of known 
Galactic RGB stars to about 40. We have obtained spectra to confirm the identifications. 
The fact that four out of the five newly identified RGB stars are "cool" (Tg// < 6000 K) 
rather than "warm" (Tg// > 6000 K) suggests that the preponderance of warm RGB 
stars among the existing sample is a selection bias. These cool RGB stars are redder 
and fainter than their warm counterparts and may have been missed in surveys done 
with blue plates. Based on the number of new RGB stars discovered in the MAGHO 
bulge fields, there may be ~250 RGB stars in the reddened "exclusion" zone toward the 
bulge. 

Subject headings: Galactic Bulge, R Goronae Borealis stars, Stellar Evolution 

1. Introduction 

R Goronae Borealis (RGB) stars are believed to be in a very short-lived evolutionary phase 
which is not yet well understood (Glayton 1996). Our ability to comprehend these stars depends 
critically on enlarging the small sample of known objects and also on understanding the severe 
selection biases at work. With irregular dramatic declines of several magnitudes over only tens 
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of days, these carbon-rich, hydrogen-deficient supergiants have proven rare, with only about 35 
RGB stars known in the Galaxy (Glayton 1996). Their spectacular declines in brightness at optical 
wavelengths are attributed to the ejection of material which condenses to form thick dust clouds 
that obscure the photosphere. 

Two competing theories for the evolution of RGB stars arc the Double Degenerate (DD) and 
the Final Helium Shell Flash (FF) models. The former involves the merger of two white dwarfs in 
a binary system (Iben, Tutukov, & Yungleson 1996a; Saio k, Jeffery 2000). In the FF scenario, a 
pre-white dwarf expands to supergiant size through a helium shell fiash (Iben et al. 1983; Renzini 
1990). This model suggests a relationship between RGB stars and PN; the evidence for such a 
relationship became stronger recently when outbursts were observed for three hot PN central stars 
(V4334 Sgr, FG Sgc, V605 Aql) that transformed them into cool giants with the spectral properties 
of an RGB star (Kerber et al. 1999; Asplund et al. 1999; Glayton k. De Marco 1997; Gonzalez et 
al. 1998). 

Recently, Alcock et al. (2001) reported the discovery of eight new RGB stars in the Large 
Magellanic Gloud (LMG) using MAGHO Project photometry. Five RGB stars have been discovered 
in the SMG (Tisserand et al. 2004; Morgan et al. 2003). Recently identified Galactic RGB stars 
include Z UMi (Benson et al. 1994), ES Aql (Glayton et al. 2002) and V2552 Oph (Hesselbach 
et al. 2003; Rao & Lambert 2003). Most of the newly discovered stars are cool (Tg// < 6000 K) 
RGB stars, implying that observational bias and not population may have precluded their discovery 
until now. Lawson & Gottrell (1990) discussed an exclusion zone towards Galactic center (GG), 
where excessive reddening prevented the discovery of Galactic bulge RGB stars in previous studies, 
such as the Harvard Observatory plate survey which had a limiting magnitude of B'~12. Alcock et 
al. (2001) suggest that My oc Tg//. So, cool, reddened RGB stars may have been systematically 
missed in previous surveys. This paper reports the results of a search of the MAGHO Project 
database for new RGB stars in the Galactic bulge. 



2. Observations 

The MAGHO Project (Alcock et al. 1992) was designed to search for gravitational microlens- 
ing events. It collected two-color photometric data on millions of stars in the Large and Small 
Magellanic Glouds, and the Galactic bulge. The survey spanned seven years, from 1992 to 1999, 
and utilized a dedicated 50-inch telescope at Mt. Stromlo, Australia. The camera was mounted at 
prime focus, had a field of view of 0.5 deg^, and was custom built for this project (Stubbs et al. 
1993). Two-color photometry in blue ( Bmacho, ~4400-5900 A) and red ( Rmacho, ~5900-7800 A) 
light was obtained simultaneously by using a dichroic beam splitter. The camera system incorpo- 
rated eight 2048x2048 pixel Loral GGDs, with 2x2 mosaics in each focal plane. The 15 /xm pixels 
mapped to 0^'63 on the sky. The photometric data were reduced using SODOPHOT, a profile- 
fitting photometry routine derived from DoPHOT (Schecter, Mateo, & Saha 1993). Galactic bulge 
data were obtained for 94 fields, covering approximately 50 deg^. The database was searched for 
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RCB star candidates as described in Section 3. MACHO can detect stars down to about 20th 
magnitude and is saturated for stars brighter than about 12th magnitude. 

We utihzed the cahbration relation of Popowski et al. (2003) (based on Alcock et al. 1999) to 
convert to Johnson V and Kron-Cousins R bandpasses: 
V = 23.70 + 0.82 Bmacho + 0.18 Rmacho 
R = 23.41 + 0.18 Bmacho + 0.82 Rmacho- 

The MACHO Hght curves, cahbrated in this way for the new RCB candidate stars are shown in 
Figure 1; {V — R)kc colors are also plotted. The stars are listed in Table 1 using the standard 
"field.tile. sequence" number, which designates a particular object in the database, and their J2000.0 
equatorial coordinates. Finder charts are shown in Figure 2. 

We obtained JHK photometry from the Two Micron All Sky Survey (2MASS) project for each 
of the stars. The data are listed in Table 2 with the Modified Julian Dates of the observations. 
The typical 1-a 2MASS errors are ± 0.03 mag in all three bands. 

Spectroscopic observations for all of the stars except 118.18666.100 were obtained in 2003 May 
from Las Campanas observatory, Chile, with the Low-Dispersion Survey Spectrograph (LDSS) 
mounted on the 6.5m Magellan II telescope. The grating used gives a resolution of ~10 A, and 
a useful range of about 3600-7000 A. The slits used were long enough to simultaneously record 
background sky and object, allowing accurate sky subtraction for each star. The spectrograph was 
wavelength calibrated using two arc lamps, and flux calibration was obtained by observing one 
standard star. Images were also acquired at Magellan in order to check the star IDs. Spectra for 
118.18666.100 were obtained on 2003 July 7 using the Steward Observatory's 2.3m Bok telescope 
at Kitt Peak, Tucson, with the B&C spectrograph. The grating gives a resolution of ~5 A, and 
a useful range of about 3200-8000 A. The spectra of the candidate stars are shown in Figure 3. 
Twelve additional candidate stars were rejected on the basis of spectra obtained at Steward and 
Las Campanas. Three of these stars belong to the AM Her class of magnetic cataclysmic variables 
(Hellier 2001). The lightcurve of an AM Her star, featuring active and inactive states, looks very 
similar to an RCB lightcurve. Few AM Her stars are identified in the optical (Drissen et al. 1994) 
so searching for RCB-like lightcurve behavior seems to be an efficient method of discovering these 
rare variables. The remainder of the rejected candidates are K and M stars. They are listed in 
Table 3. 

3. New RCB Stars in the Galactic Bulge 

The signature of an RCB star is its lightcurve (Clayton 1996). A star may stay at maximum 
brightness for months or years, then show unusual sharp declines at irregular intervals. The declines 
happen quickly, with a drop of 3 magnitudes or more taking a few days or weeks. Recovery typically 
lasts months or years (Payne-Gaposchkin & Gaposchkin 1938). The RCB stars are cool (Te// = 
5000-7000 K), hydrogen-deficient carbon stars, as revealed by their spectra. The hydrogen lines are 
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weak at best, as is the G-band (CH)^. RGB stars exhibit strong carbon Unes and molecular bands 
of GN and C2 that dominate their spectra. However, they typically have little or no ^^G. RGB 
stars often exhibit regular or semi regular pulsations with AV of a few tenths of a magnitude and 
periods of 40-100 days. They typically also have an infrared excess. 

Three previously known RGB stars lie within the MAGHO Galactic bulge fields. These are 
V739 Sgr (123.24377.5), V3795 Sgr (161.2445.2938) and VZ Sgr (117.25166.4791); however, these 
stars do not lie within the top fields, meaning there was very little data taken for them. GU Sgr 
lies just outside the area covered by the MAGHO fields. 

We searched the MAGHO database for stars with large declines in brightness, and further culled 
candidates by excluding stars whose declines were periodic. The seven year MAGHO lightcurves for 
these candidates were then examined visually to identify likely RGB stars. We used spectroscopic 
data to confirm the final candidates as RGB stars or to otherwise classify them. All of our newly- 
discovered RGB stars are present in sky surveys, such as 2MASS, but have not previously been 
identified as variable stars. We applied the criteria used by Alcock et al. (2001) to identify the 
candidates as RGB stars. The rate of decline, dm/dt (magnitudes/day), quantifies the characteristic 
sharp decline, and the lightcurves were examined by eye for evidence of pulsations. The spectra 
were examined for evidence of hydrogen by looking at the Balmer lines and the G band of GH at 
4300 A. The presence of ^^G was searched for in the isotopic bands of G2 and GN. In particular, 
the Swan bands, ^^G^^G and ^^G^^G near 4700 A, other G2 bands in the 6000-6200 A region, and 
the ^^GN band near 6250 A were examined. The photometric and spectroscopic characteristics are 
listed in Table 4. 

One star, 118.18666.100, shows similar spectral features to other warm RGB stars, such as W 
Men, although the S/N of the spectra were not sufficient for a definite identification. The typical 
low-resolution, warm RGB spectrum is almost featureless since there are no Balmer lines and the 
molecular features are very weak. The lightcurve of 118.18666.100 shows one sharp decline of > 2 
mag followed by a slow return to maximum. It experiences pulsations similar to other RGB stars. 

The four other stars, 135.27132.51, 301.45783.9, 308.38099.66, and 401.48170.2237 show spectra 
typical of cool RGB stars with very strong molecular absorption bands of G2 and GN. Two of the 
stars, 401.48170.2237 and 301.45783.9 are very active, showing several deep declines. The lightcurve 
coverage isn't as complete for 308.38099.66 and 135.27132.51 but the former shows evidence for 
two declines and the latter shows evidence for one. Pulsations can be seen in the lightcurves of 
401.48170.2237, 301.45783.9, and 308.38099.66. The fragmentary lightcurve of 135.27132.51 does 
not permit a search for pulsations. 



^V854 Ccn, one of the most active RGB stars, shows fairly strong Balmer lines and CH band (Kilkenny & Marang 
1989; Lawson & Cottrell 1989). 
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4. IR Excesses 

All of our newly-discovered RGB stars have been observed in the infrared by the 2MASS survey. 
The 2MASS data are listed in Table 2. Their positions in the H-K/J-H color-color plot would fall 
right on those observed for S Aps, a cool RGB star. See Figure 7 of Feast (1997). The 2MASS data 
only give a snapshot of each star. One star, 118.18666.100, was at or near maximum light when 
observed; 301.38099.66 was recovering from a decline; 308.38099.66 was beginning a sharp decline; 
and 401.48170.2237 was in a deep decline. The three stars in decline lie above and to the right of 
the stars at maximum light following the behavior seen for other RGB stars (Feast 1997). Due to 
the fragmentary nature of the 135.27132.51 lightcurve, we don't know where it was when observed 
by 2MASS. But from its position in the H-K/J-H color-color plot, it was not in a decline. 

5. The Population of RGB stars 

Recently, Alcock et al.(2001) estimated 3200 Galactic RGB stars, based upon the number of 
observed RGB stars in the LMC and extrapolating to the Galaxy. Both of the evolutionary theories 
referenced in the introduction, the Double Degenerate and the Final Helium Shell Flash models, 
suggest that the RGB stars arc an old population (Glayton 1996). The Double Degenerate scenario 
would suggest a population of about 1000 Galactic RGB stars (Webbink 1984). Iben, Tutokov & 
Yungleston (1996b) put the Galactic RGB population resulting from the same scenario to be about 
300 stars, and they calculate that the Final Flash scenario would imply anywhere from 30 to 2000 
RGB stars at any given time. All of the evidence thus far suggests that there are many more than 
the 40 known RGB stars in the Galaxy. 

It is a significant extrapolation from the current small sample to estimate the population of 
Galactic RGB stars. There are some severe selection biases affecting the known sample. Most 
of the RGB stars were discovered on plate surveys that went down to B~12. Discovery is also 
hampered by the fact that the RGB stars start at maximum and then get fainter so the maximum 
light brightness of the star has to be significantly above the plate limit to be detected as an RGB 
star. In the LMG where a reasonably unbiased census was possible, Alcock et al. (2001) found 
that My oc Tg// such that My falls from -5 to -3 as Tg// goes from ~7000 to 5000 K. They also 
found that the cool (Tg// < 6000 K) RGB stars far outnumber the warm (Tgjj > 6000 K) RGB 
stars. Of the Galactic RGB stars identified before 1988, primarily in plate surveys, that have Tgjj 
estimates, 16 are warm and 4 are cool. Starting with the fluke discovery of V854 Gen in 1988^, 
thirteen additional RGB stars have been identified. Of these, ten are cool and three are warm, 
just the opposite of the sample discovered earlier. These new RGB stars have a cool/warm ratio 
similar to that seen in the LMG. Therefore, cool RGB stars may have been heavily undercounted 
in previous Galactic surveys because they are redder and fainter than their warmer counterparts. 



^It is V=7.1 at mciximum light but liad spent most of tlie 20tli century in a deep decline. 
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The selection effects are exacerbated in the direction toward the GC where reddening is high. 

With the addition of the five new RGB stars reported in this work, there are now about 40 
RGB stars known in the Galaxy. Their distribution on the sky is shown in Figure 4. This larger 
sample, for which we now have more accurate estimates of Tg// and My, allow us to significantly 
improve on the distributions attempted previously (Lawson Sz Cottrell 1990; Drilling 1986). The 
new stars do not show a significantly different distribution from those discovered earlier. It is very 
hard to discern the nature of the population of the RGB stars in the Galaxy based on such a small 
sample. There is no indication of a difference in distribution for the warm and cool RGB stars. 
In Figure 4b, we have plotted the "exclusion zone" postulated by Lawson Sz Gottrell (1990) where 
RGB stars may have been missed due to large amounts of dust extinction. The addition of newly 
discovered stars plus additional stars having Tg/ / estimates have filled in the nearby portion of the 
zone but the actual distribution of the RGB stars remains mysterious. 

The distribution on the sky and radial velocities of the RGB stars tend toward those of the 
bulge population (Drilling 1986; Jeffery et al. 1987). In the discussion following Drilling's paper, 
it is reported that new data show that the radial velocities of the RGB stars are significantly lower 
than those of the EHe stars. The EHe stars are hotter than the RGB stars and don't make dust. 
The scale heights for these two groups of stars have been derived as z = 1700 pc for the EHe 
stars and 400 pc for the RGB stars (Iben & Tutukov 1985). So while the EHe stars seem to be 
bulge/Population II stars, the RGB stars may be more like old disk/Population I stars. Figure 4 
shows that while the RGB stars are concentrated toward the GG in the sky, most of them lie 2-5 
kpc from the Sun and 3-5 kpc from the GG. One can easily believe that the exclusion zone could 
be populated with very faint reddened RGB stars. But if there is a spheroidal distribution of RGB 
stars 5 kpc in radius around the GG, then there should be large numbers of lightly reddened RGB 
stars, with V mags brighter than 12, lying above and below the plane of the Galaxy. See Figure 4b. 
This is not seen. Taking into account the selection bias constraints outlined above, the distribution 
of RGB stars is consistent with a thick disk population which has a scale height of ~750 pc and a 
scalelength of 2.8 kpc (Robin et al. 1996). As shown in Figure 5, the distribution of the RGB stars 
is also similar to the distant PNe (Drilling 1986). See Figure II-l of Pottasch (1984) for comparison. 
There is a sharp peak in the density of RGB stars around the GG very much as there is in the PNe 
consistent with a Bulge population. However, the distant PNe are thought to lie close (< 1 kpc) to 
the GG and as shown in Figure 4, the known RGB stars are actually quite far from the GG, nearer 
to the Sun. 

Of the 94 MAGHO Bulge fields, only a fraction had consistently good coverage over the seven 
years. For each MAGHO field, we estimated the fraction of the mission with coverage (as opposed 
to the percentage of nights observed). Using this method, we estimated that thirty-three fields had 
coverage of less than 10% and twenty-seven fields had coverage equal to 40%. Of the five stars that 
we discovered, four lie in fields with 40% coverage. Any RGB stars brighter than V~12 will not 
be detected in the MAGHO data because they will be saturated on the GGD. This is probably a 
very small number. Alcock et al. (2001) estimated the detection rate for finding RGB stars in the 
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MACHO database to be 75%, based upon the AAVSO monitoring of 31 known RGBs. Of these, 
23 experienced a decline of greater than three magnitudes during the seven year time span of the 
MACHO project. Given that each MAGHO Bulge field contains 0.5 deg^, we are surveying ~13.5 
deg^ in the top 27 fields in MAGHO Bulge database. This is less than 3% of the area of the Lawson 
k, Gotrell (1990) exclusion zone. Based on the 5 new stars identified and assuming a 75% success 
rate, there would be ^250 RGB stars in the exclusion zone to a depth of about 7 kpc, the distance 
to the farthest RGB star discovered in this study. 

6. Summary 

Our search of the MAGHO bulge photometry database for new RGB stars had the following 
results. 

• Five new RGB stars were discovered bringing the total known in the Galaxy to ~40. 

• The discovery ratio of 4:1 for coohwarm RGB stars in this study agrees with the predominance 
of cool RGB stars found in the LMG, and indicates a selection bias for previously discovered RGB 
stars. 

• The true spatial distribution of RGB stars is still uncertain, although the projected maps indicate 
a possible thick disk population. 

• We estimate the population of RGB stars in the reddened exclusion zone toward the Galactic 
Genter to be about 250. 

To determine the true distribution and population of RGB stars in the Galaxy, MAGHO-like 
surveys will have to be done of other patches of sky at different Galactic latitudes and longitudes. 

This project was supported by the NSF/REU grant AST-0097694 and the Nantucket Maria 
Mitchell Association. Wc would also like to thank Vladimir Strelnitski, Director, Maria Mitchell 
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ter/Galifornia Institute of Technology, funded by the NASA and NSF. KHG's work was performed 
under the auspices of the U.S. Department of Energy, National Nuclear Security Administration 
by the University of Galifornia, Lawrence Livermore National Laboratory under contract No. W- 
7405-Eng-48. 
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Fig. 1. — MACHO Project lightcurves for confirmed RGB stars 301.45783.9 (upper 
left), 308.38099.66 (upper right), 401.48170.2237 (middle left), 135.27132.51 (middle right), and 
118.18666.100 (bottom). 




Fig. 2.— Finding charts for RGB stars, 301.45783.9 (upper left), 308.38099.66 (upper right), 
401.48170.2237 (middle left), 135.27132.51 (middle right), and 118.18666.100 (bottom). 
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Fig. 3.— Spectra of confirmed cool RGB stars, 301.45783.9 (upper left) ,308.38099.66 (upper right), 
401.48170.2237 (middle left), and 135.27132.51 (middle right). These spectra show strong ab- 
sorption bands of C2 and CN (Alcock et al. 2001). Spectrum of the confirmed warm RGB star 
118.18666.100 (bottom left). This spectrum shows no strong features. Spectra of typical warm (W 
Men) and cool (HV 5637) RGB stars with molecular bands marked for comparison (bottom right) 
(Alcock et al. 2001). 
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Fig. 4. — Upper Left panel: Aitoff projection showing the distribution of RGB stars on the sky. 
Open circles are new MACHO Bulge RGB stars, filled circles are RGB stars identified before 1988 
and open squares are RGB stars identified since 1988. Upper right, lower left and lower right panels 
show y vs z, X vs z and x vs y distributions of the RGB stars in units of kpc. The symbols are 
the same as in the upper left panel. The locations of the Sun and the Galactic center (GG) are 
marked. In the upper right panel, the solid lines show the exclusion zone from Lawson & Gottrell 
(1990) at b=+15° and -10° in the direction of the GG. 
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Fig. 5. — Histogram distribution of RGB stars with Galactic latitude and longitude. See text. 
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Table 1. MACHO Galactic Bulge RGB Stars 



Name^ MACHO Name a(2000) (5(2000) 



Known RGB Stars 



V739 Sgr (123.24377.5) 


MAGHO*18:13:13.1-30:15:50 


18 


13 


13.11 


-30 


15 


50.4 


V3795 Sgr (161.24445.2938) 


MAGHO*18:13:24.2-25:46:43 


18 


13 


24.23 


-25 


46 


43.3 


VZ Sgr (117.25166.4791) 


MAGHO*18:15:09.1-29:42:29 


18 


15 


09.13 


-29 


42 


28.8 


Gonfirmed New RGB Stars 
















118.18666.100 


MAGHO*17:59:52.2:-29:39:50 


17 


59 


52.23 


-29 


39 


50.0 


135.27132.51 


MAGHO*18:19:33.9-28:35:58 


18 


19 


33.87 


-28 


35 


57.8 


301.45783.9 


MAGHO*18:32:18.6-13:10:49 


18 


32 


18.60 


-13 


10 


48.9 


308.38099.66 


MAGHO*18:19:27.4-21:24:08 


18 


19 


27.36 


-21 


24 


08.2 


401.48170.2237 


MAGHO*17:57:59.0-28:18:13 


17 


57 


59.02 


-28 


18 


13.1 



'^Includes MAGHO F.T.S. number. 

Table 2. 2MASS JHK Photometry 



Name MJD-2400000 J (mag) H (mag) K (mag) 



118.18666.100 


51010 




9.2 


8.1 


7.5 










135.27132.51 


51407 




10.4 


9.7 


8.9 










301.45783.9 


51366 




12.2 


11.0 


9.9 










308.38099.66 


51339 




11.9 


10.1 


8.6 










401.48170.2237 


51010 




11.5 


9.8 


8.4 
















Table 3, 


. Final Identification 










Name 


^max 




(V-R) 


E(B-V) 


Am At 


dm/dt 


13c 


H/GH 


Pulsati 


118.18666.100 


16.6 


14.9 


1.7 


3.0 


>1.5 25 


0.06 


None 


None 


yes 


135.27132.51 


14.3 


13.1 


1.2 


1.0 


>3.5 75 


0.047 


None 


None 


? 


308.38099.66 


17.3 


15.2 


2.1 


2.8 


>5 




None 


None 


yes 


301.45783.9 


16.3 


14.8 


1.5 


1.6 


5 110 


0.05 


None 


None 


yes 


401.48170.2237 


14.5 


12.8 


1.7 


2.0 


6 100 


0.06 


None 


None 


yes 



-16- 



Table 4. Non-RCB Stars 



Name 


^max 


l^max 


(V-R) 


Probable Classification 


159.25614.43 


15.2 


13.9 


1.4 


M 


122.22949.21 


15.1 


13.7 


1.4 


M 


154.29038.66 


16.4 


16.2 


0.2 


AM Her 


149.27752.4439 


16.9 


16.7 


0.2 


AM Her 


311.37557.169 


15.3 


15.1 


0.2 


AM Her 


304.36408.33 


17.4 


15.0 


2.4 


K Star? 


305.36243.133 


16.6 


16.1 


0.5 


? 


104.20900.3659 


12.5 


13.5 


1.0 


K star 


102.26243.24 


13.6 


12.5 


1.1 


M 


303.44915.18 


15.8 


14.1 


1.7 


M 


102.23507.2965 


10.2 






M 


148.26072.30 


14.2 


13.1 


1.1 


M 


104.20385.189 


17.5 


15.5 


2 


M 



